Cooperative Behavior and Pattern Formation in Mixtures of Driven and Non-Driven 

Colloidal Assemblies 
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We simulate a disordered assembly of particles interacting through a repulsive Yukawa potential 
with a small fraction of the particles coupled to an external drive. Distortions in the arrangement 
of the nondriven particles produce a dynamically induced effective attraction between the driven 
particles, giving rise to intermittent one-dimensional stringlike structures. The velocity of a moving 
string increases with the number of driven particles in the string. We identify the average stable 
string length as a function of driving force, background particle density, and particle charge. This 
model represents a new type of collective transport system composed of interacting particles moving 
through deformable disorder. 

PACS numbers: 82.70.Dd 
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The collective transport of interacting particles driven 
through disordered backgrounds has been studied exten- 
sively in systems such as moving vortex lattices in the 
presence of random pinning .1), driven charge density 
waves , and sliding friction p| ■ A variety of nonequi- 
librium behaviors occur in these systems, including frac- 
tal flow patterns, avalanches, nonlinear velocity- force re- 
sponses, and dynamic reordering transitions. The collec- 
tive effects arise due to the particle-particle interactions, 
while the disorder in the substrate is quenched and does 
not change with time. It is also possible to have sub- 
strates which can distort and respond to the motion of 
the driven particles, such as in the case of a single col- 
loidal particle driven through a disordered background 
of other particles which are not coupled to the external 
drive Both simulations [f| and experiments Q 

show that such a system can exhibit a nonzero threshold 
force for motion as well as a nonlinear velocity vs applied 
force response. The driven particle moves in an intermit- 
tent manner and causes local rearrangements and distor- 
tions in the surrounding bath of colloids. An open ques- 
tion is what happens if there are multiple particles driven 
through a background of nondriven particles, rather than 
only a single driven particle. In particular, it is possible 
that the ability of the surrounding media to react to the 
motion of driven particles could induce an emergent ef- 
fective driven particle-particle interaction, leading to new 
types of moving structures or to pattern formation that 
is not observed in driven systems with fixed or quenched 
background disorder. 

Recent studies have focused on a related system in 
which two sets of identical repulsive particles are driven 
in opposite directions 0, 0, El El- Here, a laning phe- 
nomenon occurs, and the species segregate into multi- 
ple streams flowing past one another. The oppositely 
driven species model has been studied both in terms of 
pedestrian dynamics, where the particles represent peo- 
ple moving in opposite directions as well as in the 



context of colloidal or charged Yukawa particles in the 
presence of a driving field [jj. The colloidal systems are 
particularly attractive for further study since binary mix- 
tures of oppositely charged colloids which can undergo 
laning transitions in the presence of a driving field have 
been produced recently in experiment This type of 
system opens a wealth of new experimental possibilities. 

To our knowledge, the case of a small fraction of par- 
ticles driven through a responsive but undriven back- 
ground has not been studied previously. In this work, 
we show that when multiple particles move through a 
background of nondriven particles, an effective emergent 
attraction arises between the driven particles despite the 
fact that all of the pairwise interactions between the par- 
ticles are repulsive. We find that intermittent moving 
strings form which have an average stable length that is 
a function of the driving force, system density, and parti- 
cle charge. Strings that are longer than average are very 
short lived. The string velocity increases monotonically 
with the string length. 

We consider a two-dimensional system with periodic 
boundary conditions in the x and y directions contain- 
ing N c particles interacting via a screened Coulomb or 
Yukawa interaction. Only a small number Njj of the 
particles couple to an external drive. The particles move 
in an overdamped background, and we neglect hydrody- 
namic effects, which is reasonable for the low volume frac- 
tion limit considered here. We use molecular dynamics 
(MD) to numerically integrate the overdamped equations 
of motion for the particles, given for particle i by dxi jdt = 
Fj C + Ff -(- F T . The colloid-colloid interaction force term 
F^ c = — qi ^iV{ r ij) contains the Yukawa interac- 
tion potential V^ry) = {qj/\ti — r$\) cxp(— «|r{ — r 3 |), 
where r^y) is the position of particle q^j) is the 

charge of colloid and all colloids have the same sign 
of charge so that all interaction forces are repulsive. The 
screening length 1/k is set to 1/2 in all our simulations. 
The background colloids are composed of a 50:50 mix- 
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FIG. 1: Images of the system composed of a disordered array 
of background colloids (small circles) and particles that couple 
to the external drive (large circles) for qo /q = 6.67, n c = 0.94, 
F D = 3.0, and N D /N C = 0.01. (a) Initial configuration; (b) 
the same system after 10 5 MD steps; (c) a typical steady state 
snapshot after 10 MD steps; (d) snapshot of a string breaking 
into shorter strings. The driving force is from left to right. 

ture of two charges with qi/q2 = 1/2, which produces 
a noncrystalline background arrangement. The average 
background charge q = (qi + <?2)/2. For the non-driven 
particles, the constant drive force term F^, = 0. The 
colloids that couple to the external field have charge qn 
and experience a constant drive Fjj = Fox applied af- 
ter the system is equilibrated by simulated annealing. 
The thermal force Ff has the properties (Ff ) = and 
(Ff (t)F](t')) = 2Trr]kBT8i j S{t - t 1 ). For a system of 
length L the colloid density is n c = N c /L. We consider 
both the constant n c case, where L and N c are increased 
simultaneously, as well as the case where n c is increased 
by fixing L and increasing N c . 

In Fig. 1 we show four snapshots which highlight 
the dynamically induced attraction of the driven par- 
ticles, indicated as large black dots, and the formation 
of the string structures. Here the density n c = 0.94, 
QD/q = 6.67, and Fd — 3.0. Fig. 1(a) shows the ini- 
tial state in which the driven particles are well separated 
from each other. In Fig. 1(b), after 10 5 MD time steps, 
the driven particles begin to form string like structures 
of length N s = 2 to 3 colloids which are aligned in the 
direction of the driving force (left to right in the figure) . 
In this noncquilibrium system, despite the fact that all 
of the pair-wise interactions are repulsive, an effective at- 
traction emerges between the driven particles. What is 
difficult to convey through still images is the fact that 




FIG. 2: (a) Time series V x (t) for a single driven colloid from 
the system in Fig. 1(a) with Nd = 20. (b) The corresponding 
time series for V y . (c) Solid line: The normalized histogram 
of V x from the complete time series for the system in (a). 
Dotted line: Histogram of V x for the case Nd = 1- (d) Solid 
line: The normalized histogram of V y for the total time series 
from the system in (b). Dotted line: Histogram of V y for the 
case Nd = 1. 

as N s , the number of particles in a string, increases, 
the average string velocity (V^f) increases so that strings 
move considerably faster than isolated driven particles. 
In Fig. 1(c) after 10 6 MD time steps, the system has 
reached a noncquilibrium steady state. Here the average 
stable length of a string (N s ) = 5 to 6. Strings with 
N s > 6 form occasionally but generally have a very short 
lifetime. In Fig. 1(c), strings with N s = 6, 5, 3, and 1 
are present. All of the strings have an intermittent char- 
acter, in that driven particles in the string can become 
detached from the string, and new driven particles can 
join the string. Typically strings add or shed one driven 
particle at a time; however, long strings with N s = 7 or 
8 tend to break up into separate strings with N s = 3, 4, 
or 5, as illustrated in Fig. 1(d). 

We note that if Nd/N c is increased while all other pa- 
rameters are held fixed, then {N s } remains unchanged 
at (N s ) = 5 to 6. We find the formation of stable 
strings over an extensive range of parameters, including 
all qo/q> 1-0 (where we have checked up to qo/q = 40), 
for n c > 0.17, and for arbitrarily large Fd- The value 
of (N s ) depends on the specific choice of parameters; 
however, the phenomenon of string formation is very ro- 
bust provided that the driven particles can distort the 
arrangement of the background particles. In Fig. 1, 
Nd/N « 0.01, so that the initial average distance be- 
tween driven particles is many times larger than the 
pair-interaction length. Thus, the attractive force be- 
tween driven particles is mediated by the distortion field 
produced in the background particles. If the fraction 
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FIG. 3: (V^f) vs the number of particles in the string N s . 
(a) The same system as in Fig. 1(a) with qo/q = 6.67 for 
Fd = 3.0 (squares), 4.0 (triangles), and 5.0 (circles), (b) The 
same system as in (a) for qo/q = 1-33 and Fd = 0.5 (squares), 
0.75 (triangles), and 1.0 (circles). 

is increased as much as an order of magnitude up to 
Nd/N c — 0.1, the preferred length of (A^) = 5 to 6 per- 
sists. We note that if /N c is too high, finite size effects 
due to the boundary conditions become important. 

In order to illustrate the intermittent behavior and the 
velocity increase as the number of colloids in a string 
increases, in Fig. 2(a) we plot a typical time series of the 
velocity in the direction of drive V x for one of the driven 
colloids from the system in Fig. 1(a). Here V x (t) shows 
a series of well defined increasing or decreasing jumps, 
with a roughly constant velocity between the jumps. The 
lower value that V x takes is around 0.39 corresponding 
to the average velocity of a driven particle that is not in 
a string. There arc several places where V x — 1.45 which 
corresponds to the average velocity of a string with N s = 
5. Additionally, there are some plateaus near a value of 
V x = 0.7 which corresponds to A^ s = 3. Fig. 2(a) also 
shows that for short times, V x > 1.5, which corresponds 
to strings with N s > 6. This figure shows that the driven 
colloid is exiting and joining strings of different lengths. 

Fig. 2(a) is only a portion of a much larger time se- 
ries which covers 10 7 MD steps. In Fig. 2(c) we plot the 
histogram of this entire time series along with the his- 
togram of V x (dotted line) for a case where there is only 
one driven particle in the system so that strings do not 
form. The N u = 1 curve shows a single peak in P{V X ) 
near V x = 0.35, which falls near but slightly below the 
first peak in P(V X ) for the multiple driven particle sys- 
tem. Thus, the lowest value of V x in Fig. 2(a) corresponds 
to time periods when the driven colloid is moving indi- 
vidually and is not part of a string. Additionally, there 
is a small peak in P(V X ) for the Nd — 20 system near 
V x = 0.8, which corresponds to the formation of strings 
with A^ s = 2. There is a large peak in P{V X ) around 
V x = 1.45 corresponding to A^ s = 5 and 6. For V x > 1.5 
there is a rapid decrease in P(V X ). This result shows 
that there is a preferred string length of (N s ) « 5 for the 
Nd = 20 system. The average velocity of a driven parti- 
cle in the N u = 20 case is much larger than the average 
velocity in the Nd — 1 case, even through the fraction of 
driven particles is only Nd/N c = 0.01. 




FIG. 4: (a) The most probable number of particles in a 
string (N s ) vs Fd for the system in Fig 1. (b) (N s ) vs qo for 
the same system as in (a) for Fd = 1-0. (c) (N a ) vs n c for 
qo = 6.67 and Fd = 2.0. (d) (V^f) vs system size L for the 
system in Fig. 1 for Fd = 1.5 (triangles), 3.0 (squares), 4.5 
(diamonds), and 6.0 (circles). 

Driven particles traveling in strings show a more pro- 
nounced transverse wandering motion than individual 
driven particles. In Fig. 2(b), we plot the time series for 
the velocity transverse to the driving force, V y , for the 
same particle as in Fig. 2(a). Here (V^) = 0, and there 
are numerous spikelike events in the positive and nega- 
tive directions. In Fig. 2(d) we plot the histogram of the 
entire time series of V v (solid line) along with the same 
histogram for a system containing only a single driven 
particle (dashed line). The Nd = 1 curve shows a well 
defined Gaussian distribution of P{V y ) centered at zero. 
For Nd = 20, the spread in V y is much larger and the 
fluctuations are non-Gaussian, as indicated by the pres- 
ence of large tails in the distribution. There are an excess 
of events at larger values of \V y | which result from the fact 
that the particles in a chain move in a correlated fashion. 
The non-Gaussian statistics also indicate that the strings 
have transverse superdiffusive behavior. 

We next examine the average velocity of the strings 
(Vjf) vs N Sl which we plot in Fig. 3(a) for systems with 
the same parameters as in Fig. 1 at Fd — 3.0, 4.0, and 
5.0. The bottom curve, which corresponds directly to 
the system in Figs. 1 and 2, shows that for N s = 1.0, 
(V x s ) = 0.35, while for N s = 5 and 6, (V x s ) w 0.15. 
The string velocity increases rapidly with N s for N s < 5 
and then increases much more slowly with N s . The 
longer strings become considerably more winding and 
experience additional frictional drag, whereas the shorter 
strings move rigidly. Due to the meandering of the longer 
strings, there is a tendency for the driven particles at the 
back of the string to be ripped away from the string very 
quickly. For curves with higher Fd , there is still a mono- 
tonic increase in the average velocity with N s . In order to 
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show that this phenomenon is robust over a wide range of 
parameters, in Fig. 3(b) we plot (V^f) vs N s for a system 
with q D /q = 1.33 for F D = 0.5, 0.75, and 1.0. For these 
parameters we again observe the string formation and an 
increase in (V^f) with N s . The best functional fit we find 
is (Vx) x m (Af s ); however, other forms can also be fit, 
such as two linear fits for the system with qu — 1.33. 

In Fig. 4(a) we show the most probable number of par- 
ticles (N s ) in a string for the same system in Fig. 1 for 
varied Fd- Here (N s ) goes through a maximum as a 
function of Fd- As Fd decreases below F& < 0.2, the 
distortion created by the driven particles in the surround- 
ing background particles decreases, and therefore there is 
less attraction between the driven particles. At high Fd, 
the velocity fluctuations along the string become impor- 
tant, which tends to limit the string length. In Fig. 4(b) 
we show that (N s ) for a system with a fixed Fd = 1.0 also 
goes through a maximum as a function of qo /q- For small 
qo/q, the driven particles produce very little distortion 
in the surrounding medium so there is little attraction 
between them. For large qo/q, the repulsion between the 
driven particles becomes more prominent and the aver- 
age velocity drops, which again reduces the distortion 
field. In Fig. 4(c) we plot (N s ) vs the particle density n c 
for a system with qp = 6.67 and Fjj — 2.0. At low n c , 
all the particles are far apart and the distortion in the 
bath particles is reduced or absent. As n c increases, the 
particle-particle interactions arc enhanced and the string 
can grow longer. Presumably for very high n c , {N s } will 
decrease again as the average velocity will have to be re- 
duced. We are not able to access particle densities in this 
range. 

In Fig. 4(d) we plot (Vjf) vs the system size L for 
Fd = 1.5 (lower curve), 3, 4.5, and 6.0 (upper curve). In 
each case, (V^) initially decreases with L but reaches a 
saturation value by L « 40. The change in (V*) at small 
L occurs because, in small systems, the periodic bound- 
ary conditions allow driven colloids to interact strongly 
with their own distortion trails. For larger Fd, larger sys- 
tem sizes must be used to avoid this effect. In addition 
to the increase in (V^) in the small systems, (N s ) also 
increases for the smallest systems. This emphasizes the 
fact that sufficiently large systems must be considered in 
order to avoid finite size effects. We have used L — 48 
throughout this work, which is large enough to avoid fi- 
nite size effects for the quantities we have measured. 

The attraction between the driven particles is medi- 
ated by the distortion in the surrounding bath of col- 
loids. The concept of noncquilibrium depletion forces, 
which are highly anisotropic and have both an attractive 
and a repulsive portion, was introduced recently for a 
case where two large colloidal particles have an additional 
bath of smaller particles flowing past them 10] . The sur- 
rounding particle density is higher in front of the large 
particle and lower in back, which causes an anisotropic 
attraction force between the large colloids that aligns 



them with the direction of flow. Although this model 
was only studied for two particles, it seems to capture 
several of the features that we find in our simulations, 
including the alignment of the strings in the direction of 
drive. This suggests that noncquilibrium depletion force 
models could be applied for multiple large particles. The 
longer strings move faster since the front particle has ad- 
ditional particles pushing on it from the back. Since the 
chains are aligned, the cross section with the background 
particles remains the same as for a single driven particle; 
however, the string is coupled N s times more strongly to 
the external field. As longer strings become increasingly 
meandering, the cross section with the background par- 
ticles increases and limits the speed that long strings can 
attain. 

In conclusion, we have investigated a system of repul- 
sively interacting Yukawa particles which form a disor- 
dered background for a small fraction of Yukawa particles 
that couple to an external drive. In this system, despite 
the fact that all of the pair-wise interactions are repul- 
sive, there is an emergent effective attraction between the 
driven particles which results in the formation of one- 
dimensional intermittent stringlike structures. The ve- 
locity of a string increases with the number of particles 
in the string. The attraction between the driven particles 
is mediated by the distortion of the surrounding particles 
and is likely a result of the recently proposed nonequilib- 
rium depletion forces. Longer strings move faster than 
shorter strings due to the additional pushing force from 
the back particles in the string. This system represents a 
new class of collective transport where multiple interact- 
ing colloids are driven through a deformable disordered 
background, which is distinct from the case where the 
background disorder is quenched. 
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